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ABSTRACT 



> 

Motivated by a recent study by Lazarian and Draine, which showed that 
a high degree of grain alignment of the paramagnetic dust is achievable if the 
rates of internal relaxation are controlled by the Barnett relaxation process, 
we undertake a study of an alternative mechanism of internal dissipation, 
namely, the inelastic dissipation of energy in oblate dust grains. We find 
that deformations at double frequency that were disregarded in earlier studies 
dominate the inelastic relaxation. Our results indicate that for sufficiently 
prolate and oblate grains, e.g. for grains with 4:1 axis ratio, or/and grains 
formed by agglomeration inelastic relaxation dominates the Barnett relaxation 
within large (a > 0.1 fim) grains. For grains with axis ratio less than 1:2 the 
^ ■ inelastic relaxation is dominant for suprathermally rotating grains. 

Subject headings: Magnetic fields; Galaxies: magnetic fields, ISM: molecular 
clouds, magnetic fields 



1. Introduction 

The problem of inelastic dissipation of a freely rotating body emerges naturally in the 
context of interstellar grain alignment. 

It is well known that the polarization of starlight is caused by the alignment of 
interstellar dust grains. This was discovered by Hiltner and Hall (see Hall 1949 and 
references therein) who showed that the polarization was entailed by the linear dichroism, 
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i.e., by the differential extinction, by clouds of aligned grains, of electromagnetic waves of 
different polarizations. 

Various mechanisms of grain alignment (see table 1 in Lazarian, Goodman & Myers 
1997) depend on the internal dissipation of energy (Purcell 1979). This dissipation aligned 
the axis of the maximal moment of inertia and the angular momentum. This changes 
grain dynamics and should be accounted for in the theory of grain alignment. Therefore 
the grain alignment depends on the internal relaxation. In some cases, e.g. the alignment 
of thermally rotating grains, the alignment is not sensitive to the exact rate of internal 
relaxation. However for grains rotating suprathermally, more detailed knowledge of the 
internal relaxation is usually required (see more in the discussion section of this paper). 

It is normally believed that the Barnett dissipation is the dominant process of internal 
relaxation within interstellar grains (Purcell 1979, Spitzer & McGlynn 1979, Lazarian 1994, 
Draine 1996). It happens due to the time-dependent (due to the precession of angular 
velocity in grain body axes) remagnetisation of the sample by means of the Barnett effect. 
Indeed, the Barnett effect develops magnetization parallel to the angular velocity of a 
rotating body because a share of the angular momentum gets transferred to the spins 
(see Landau and Lifshitz 1984). In the cause of the angular- velocity precession about the 
angular momentum (in accordance with the Euler equations), the magnetization direction 
is also precessing in the body-axes frame, and this entails dissipation. A recent study of 
paramagnetic alignment of suprathermally rotating grains (or "Purcell alignment") has 
shown that grains can be nearly perfectly aligned by the interstellar magnetic field if the 
internal dissipation is controlled by the Barnett effect (Lazarian & Draine 1997). 

Another process responsible for dissipation, known as inelastic relaxation happens as 
time-dependent stresses emerge in a body when it rotates about an arbitrary axis. In this 
paper we subject this process to scrutiny, as the existing treatment (Purcell 1979) has a 
number of difficiencies that we discuss in the paper. 

A rigorous treatment of internal relaxation is important: if the inelastic relaxation 
is more efficient than the Barnett relaxation then it is the former, not the latter, that 
determines the dynamics of crossovers and, thereby, the degree of alignment achievable by 
suprathermally rotating grains (see Lazarian & Draine 1997). 

In what follows we briefly discuss the rotation of an oblate grain (Section II); calculate 
the acceleration experienced by a point inside the grain (Section III). Then we compute the 
stresses caused by grain precession (Section IV), the rate of energy dissipation (Section V) 
and the rate of internal-dissipation-caused alignment of the major-inertia axis of the grain 
to its angular momentum. We also discuss whether the obtained formulae could be easily 
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altered in the case of prolate grains (Section VI). Our results are discussed in Section VII. 

2. Notations and Assumptions 

We consider a freely rotating cosmic-dust grain, employing two Cartesian coordinate 
systems, each with its origin at the center of mass of the body. The inertial frame (X, Y, 
Z), with unit vectors ex, ey, ez, is chosen so that its Z axis is parallel to the (conserved) 
angular momentum J (see Fig. 1). Coordinates with respect to this frame are denoted by 
the same capital letters: X, Y, and Z. We also use the system which we call body frame 
and associate with the three principal axes of inertia: 1, 2, and 3, with coordinates x, y, 
and z along these. The appropriate unit vectors are ei, e 2 , e 3 . 

Without loss of generality, one may take I 3 > Ii and 7 3 > I 2 , where Jj are the principal 
moments of inertia. 

Following Purcell (1979), we denote the angular velocity by O and reserve u for the 
rate of precession. It is a well-known fact that, as a perfectly elastic body rotates, the 
extremity of the angular velocity vector O describes a curve that is the intersection of two 
ellipsoids. One of those comes from the angular-momentum conservation: 

if VL\ + If tt 2 2 + II ttj = J 2 = const , (1) 

where ^1,2,3 are the body-frame-related components of Q). Another one, known as the 
Poinsot ellipsoid, is defined by the kinetic-energy conservation: 

httj + I 2 tt 2 2 + I 3 tt 2 3 = J Q = 2 E rot = const . (2) 

In the presence of inner dissipation, the Poinsot ellipsoid will be distorting in the cause of 

— 1/2 

time, and will eventually degenerate into a line of length I 3 oriented along the axis of 
maximal moment of inertia. That will happen after the body aligned its angular velocity 
along its maximal-inertia axis: Qi = Q 2 = 0, |f2 3 | = |J|// 3 - This mode minimizes the 
kinetic energy, with the angular momentum J fixed. 

In what follows we model the grain by an oblate symmetric top, e.g. a disk 

h > h = h = I ■ (3) 

In the grain-alignment studies such an approximation is customary (Purcell 1979, Spitzer 
& McGlynn 1979). For such a body the Euler equations read: 

sii + ^-p^ n 2 n 3 = , (4) 
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n 2 + fi 3 fii = o , (5) 

ft 3 = , (6) 

where is the body-frame components of the angular velocity. As well known, the angular 
velocity f2 nutates about the principal axis 3 at a constant rate 

u = (h - l)fi 3 , h = h/ 1 . (7) 

However from the point of view of the inertial observer it is rather the principal axis 3 
that is precessing about J (which is conserved in the inertial frame). Therefore angle 9 is 
constant. Now, let a be the angle described by the precession of vector ft about axis 3 (or 
of axis 3 about O). It is easy to see that 

tan# = (I |S1| sin a)/(I$ |f2| cos a) = h tana ; (8) 

hence a is also constant. As Q3 = Q cos a, the angular velocity of precession given by ([?]), 
except a special case of ft and J orthogonal or almost orthogonal to the maximal-inertia 
axis 3, is typically of the order of |f2|. Hence one may call the rotation and precession 
"fast motions" , implying that the relaxation is a slow process (which means that the rate 
of alignment: S<w). 

Before pursuing with the calculational part of our article, let us once again dwell 
upon the physical mechanism that changes Q3 over the spans of time comparable with the 
internal relaxation time. In terms of energy, everything is clear: the inelastic body aligns 
so as to minimize its rotational energy and still to conserve its angular momentum. This 
obvious mechanism does not seem to be instilled in the Euler equations (£|), @, @, for the 
torques remain to be zero, and nothing seems to be able to shift the solutions ^1,2,3^) from 
the form that yields a constant. In reality, however, this would-be contradiction is easily 
resolved if one recalls that inelastic relaxation stems from the inner displacements of grain 
material. These displacements result in changes of grain moments of inertia. Therefore 
in the cause of the grain's wobbling, its principle moments of inertia fluctuate, and terms 
flili 7^ should be accounted for in Euler equations. Inelasticity results in a phase shift 
between the angular velocity and /j and this causes internal alignment^. 

Since Ii depend on the grain elasticity, while the phase offset of Jj and Q is determined 
by inelastic effects, then the parameters that determine elasticity and dissipation within the 
material should enter our final formulae for the inelastic relaxation. 



1 The latter statement is equivalent to the statement that there is a phase shift between the stresses that 
act upon the grain, and the deformation of the grain. 
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3. Acceleration of a Point inside an Oblate Grain ( I 3 > Ij = I 2 ) 

Our goal is to calculate 9, the rate of the maximal- inertia axis' approaching the 
direction of angular momentum J. To achieve this goal, one has to know the rate of energy 
losses entailed by the inelastic deformation. To calculate the deformation, we shall have to 
know the acceleration experienced by a particle located inside the grain at a point (x, y, z). 
Note that we address the proper acceleration, i.e. that with respect to the inertial frame 
(X,Y,Z), but we shall express it in terms of coordinates x, y and z of the body frame 
(1,2,3). The fast processes (revolution and precession of a symmetric oblate grain) are 
described by the Euler equations (f| - §P whose solution, in neglect of the slow relaxation, 
will read 

Qi = Qj_cosut , Q 2 = fij_sinu;i; , SI3 = const (9) 

where 

Q± = Q sin a , f2 3 = Q cos a (10) 

and 

^_i_/^3 = t an o. — h tan 9 . (11) 

Besides, we shall need formulae connecting the components of fl with the absolute 
values of the angular momentum: 

tt 3 = ^ = ^ cos 9 , Ol = ^- h sin 9 (12) 

We denote the position, velocity and acceleration relative to the body frame as: r, v, a , 
while those related to the body frame (1, 2, 3) will be called r", v", a" , where we shall 
keep in mind that r = r". The acceleration a in the inertial frame looks: 

a = a" + Q, x r" + 2 fl x v" + fl x (fi x r") , (13) 

where the first and the third terms vanish, as the values of relative deformations 51/1 are 
minute in a solid body. The values of v" ~ 51 /r and a" pa 51 /r 2 , where r is the period of 
rotation, are negligible compared to velocities and accelerations in the inertial system of 
reference that by the order of magnitude are, respectively l/r and 1/t 2 . The second and the 
fourth term in Eq. ([13]), referred to the grain frame with unit vectors (ei,e2, 63), will read: 

O x r" = ei (O2 z — (I3 y) + e 2 ((I3 x — £l\ z) + e 3 (tti y — fi 2 x) 

= ei uj z Qi + e 2 oj z VL 2 + e 3 ( — u y Q 2 — w % ^1) , (14) 

and 



n x (n x r ") = n (n • r") - r" n 2 = 
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ei {fix ( tt 1 x + Q 2 y + ^3 z ) - x Q 2 } + 
e 2 {Q 2 (Q 1 x + tt 2 y+ Q 3 z) - ytt 2 } + 
e 3 {fi 3 ( ^ x + n 2 y + Vl 3 z ) - zQ 2 } . (15) 

All in all, 

a = ei { — x Vt 2 2 — x Vt 3 2 + y Vt 2 Vti + z Vt 3 Vti + z to Q± cos cut } 
+ e 2 { x £l 2 Cli — y fix 2 — y fi 3 2 + z Sl 3 f2 2 + Q± sin cut } 
+ e 3 { x £li f2 3 + y f2 3 f2 2 — -2 fij_ 2 — fi^ (x cos cut + y sin cut) } . (16) 

We had to present our calculation in great detail because the latter formula considerably 
differs from the one in (Purcell 1979), eq. 28. To be doubly sure of our formula being correct, 
we derived it also in the intermediate, "Eulerian", system of reference (see Appendix |A|). 



4. Stresses and Strains Caused by the Precession 

In this section we derive the stresses and strains produced in the rotating body by the 



time-dependent terms in (|T6), as only these terms influence the rate of energy dissipation 
in the grain. 



With aid of (^) one can easily split fll6|) into a time-independent and time-dependent 
parts: 

a = a + a t , (17) 



where 



and 



a = - ( ex x + e 2 y ) ( Q 2 3 + ^Q 2 ± ) - e 3 z Q] 



a t = e i { 2 ^J- x cos 2 ^ ^ S ^ n Z ^ 3 ^ C ° S ^ ^ 

+ e 2 {— Q 2 ± x sin 2ut — — fi^ y cos 2u;t + 2 Q± Q3 h sin cut } 

+ e 3 {£l± fi 3 (2 — /i) (x cos Lut + y sin ujt ) } . (19) 

Now we see that dissipation of energy will be taking place at two modes one of which will be 
of double frequency. Hereafter we shall be taking into consideration only time- dependent, 
i.e. produced by a t , inputs in the stresses and strains. 
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Following Purcell, we shall model the oblate grain by a prism of sizes 2a x 2a x 2c, 
(c < a). The stresses vanishing on the boundaries of our rectangular prism readQ: 



oQ 2 

!-^(x 2 - a 2 ) cos2ujt 



a 



yy 



(y 2 — a 2 ) cos2ivt 







>xy 



a, 



a 



j n 2 ± (x 2 + y 2 - 2a 2 ) sin 2ut , 
h (z 2 - c 2 ) + (2 - h) (x 2 - a 2 
h (z 2 - c 2 ) + (2-/i) (y 2 - a 2 



cos cut 



sin cut 



which leads to forces per unit volume: 

Fx di(J x i 

p 



^Q± 2 x cos2u;t + ^Q± 2 y sin2c<jt 



pQ,±_Q,$hz cos cut 



x sin 2^ — 2^- L y cos ^ UJ ^ + pQ±Q-zhz sin tut 



(20) 

(21) 
(22) 
(23) 

(24) 

(25) 
(26) 



F z = diO Z i = pf2^f2 3 (2 — h)x cos cut + pfi^fi 3 (2 — h)y sin cut , 

in full accordance with (|19"D. Here p is the density of the grain material. For a 2a x 2a x 2c 
prism, the moment of inertia i3 and the parameter h are function of the half-sizes a and c: 



16 



p a c 



and 



* = 7 



1 + (c/a) 2 
According to (3.12), (3.13), 



h 



1 + (c/a) 2 



fi 3 



fin COS0 



fi. 



fin /i sin# 



(27) 
(28) 

(29) 



fi = J// 3 being the typical angular velocity of a grain. Our knowledge of the stresses 
should be sufficient for computing the rate of energy losses in the body. 

Our expressions for the stress-tensor components differ very significantly from the 
appropriate expressions derived in (Purcell 1979). There are two reasons for it. The 
first reason is the afore-mentioned Purcell's miscalculation of the acceleration experienced 
by a point inside the wobbling body. (Because of that Purcell completely missed the 



2 Our formulae for the stresses differ from those in Purcell (1979). Purcell's expressions do not include 
the double-frequency terms, and they do not vanish on the boundaries of the grain, which they should do. 
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double-frequency contribution to the stress tensor. Later we shall show that this input is 
of the leading order.) The second reason is that Purcell forgot to impose proper boundary 
conditions upon the stress-tensor components: his stresses fail to vanish on the boundaries 
of the body. 

Our goal now is to calculate the average deformation- caused energy, stored in the 
tumbling body, and to estimate the energy dissipation, using an appropriate quality factor. 
We shall take into account that the deformation of the grain is neither purely elastic 
nor purely plastic, but is a superposition of the former and the latter. It is then to be 
described by the tensor e# of viscoelastic strains and by the velocity tensor consisting of 
the time-derivatives e^. The stress tensor will now be separated into two components: the 
elastic stress and the plastic (viscous) stress: 

a ¥ = af +af , (30) 

where the components of the elastic stress tensor are interconnected with those of the strain 
tensor (Landau and Lifshitz 1976): 

6ij = 6u T L^fL + ( a f _ I 5 Tr a&) — , (31) 



a\f = KS^ Tre + 2 // U, ~ ^ S tj Tre) , (32) 



1 

3 

\i and K being the isothermal shear and bulk moduli, and Tr standing for the trace of a 
tensor. Components of the plastic stress are connected with the strain derivatives as 

% = ^ + - I * Tr aW) i , (33) 



» 1 1 



<jf = CSij Tre + 2r, fe y - - 5 {j Tre) , (34) 
where rj and ( are the shear and stretch viscosities. 

5. Dynamics of a tumbling body 

The kinetic energy of a rotating body reads, according to @, ©, @, and fll2]) : 



E rot = x - (h n? + h n 2 2 + h fig) = \ I 1 ^1 + h fi 3 2 



- sin 2 9 + -J- cos 2 9 
1 h 



J (35) 
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wherefrom 



dE rn t J 2 



(h — 1) sin 6* cos 9 = lu J sin 9 . (36) 

Formula (j3~6| ) provides an insight in how the change of rotational energy yields a change of 
9. If we calculate the rate of energy losses, E rot , it will be easy to find the rate of alignment 
9 as (dE rot /d9)~ 1 E rot . How to compute E rot l The rotational energy changes via the 
inelastic dissipation, so that 

E rot = W (37) 

Hence what we have to find is the rate of the elastic-energy losses W. Then, with aid of 
(|36|) we shall calculate the rate of alignment: 

d6 ( dE rot \ dE rot , . ,_i • , , 

it = {-ir) ~ir = {ujJ sm6) w m 

In our case, dissipation is taking place at two modes: 

W = WW + W'^) = lu — - + 2uj — — w - \W {UJ) + 2W {2uj) \ (39) 
where we used the fact that the quality factor is almost frequency-independent: 

Q(2U,) w Q(U>) w Q _ 



6. Rate of Energy Dissipation at low temperatures 

Before pursuing to our calculations of the energy-dissipation rate W , several prefatory 
notes will be in order. As well known, at low temperatures materials are fragile: when the 
deformations exceed some critical threshold, the body will rather break than flow. At the 
same time, at these temperatures the materials are elastic, provided the deformations are 
beneath the said threshold: the sound absorption, for example, is almost exclusively due 
to the thermal conductivity rather than to the viscosity. These facts may be summarized 
like this: at low temperatures, the viscosity coefficient 7] has, effectively, two values: one 
value - for small deformations (and this value is almost exactly zero); another value - for 
larger-than-threshold deformations (and that value is high^). 

At high temperatures materials become plastic, which means that the shear viscosity 
7] gets its single value, deformation-independent in the first approximation. On the one 



3 Effcctively it may be put infinity because, as explained above, the body will rather crack than 
demonstrate fluidity. 
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hand, this value will be far from zero (so that the scattering of vibrations will now be 
predominantly due to the viscosity terms, not due to the thermal conductivity). On the 
other hand, this value will not be that high: a plastic body will rather yield than break. All 
this is certainly valid for the stretch viscosity ( as well. 

Since at low temperatures the bodies manifest, for small displacements, no viscosity 
( wq ~ u( <C n ~ K ), the stress tensor will be approximated, to a high accuracy, by its 
elastic part: instead of the system ([30]) - ([34]) we shall simply write: 

Tr a f 1 \ 1 

€lJ = Sij ~9~K + ~~ 3 Sij Tra ) Yn ' ' !U) 

This will enable us to derive an expression for the elastic energy stored in a unit volume of 
the precessing body: 

1 1 J/2/i 1 



dW/dV = £ cy ay = - A \ ( - \ ) (Tr a) 2 + a l3 a {j 



{ - \ (Tr a) 2 + a 2 xx + a 2 yy + a 2 zz + 2 (a% + a 2 yz + a 2 2 ,)} (41) 



i r i 

4^ 



where we made use of the expressions connecting the shear and bulk moduli with the 
Young modulus E and Poisson's ratio a: since K = E/[3(l — 2a] and fi = E/[2(l + a] 
then 2fi/(9K) — 1/3 = —a/(l + a). As for frozen solids Poisson's ratio a is typically 
about 0.25, we put 2fi/(9K) — 1/3 ~ — 1/5 . Mind that the body is trembling at two 
frequencies: ui and 2 uj . Anticipating the different rates of dissipation at these two modes, 
we shall split the total elastic energy into two parts: 

dW/dV = dW^/dV + dW {2ul) /dV (42) 

where, according to (pfi| ) - (p3|), 

dW^/dV = ± {a 2 yz + a 2 zx } (43) 

dW^/dV = i- {-\{Traf + a 2 xx + a 2 y + 2^} . (44) 

^From now on we shall be interested in the energies averaged over several periods of the 
precession. Therefore we shall substitute sin 2 ... and cos 2 ... by 1/2, and shall omit 
expressions sin . . . cos . . . , sin uot sin 2ut and cos ut cos 2ut . With the above reservation 



and 



being beared in mind, expressions fl20|) - (|23|) for the stresses will yield: 



(Tra) 2 = [ Pj f] ^ { X \ - ytf \ , (45) 
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at = I^Y «« fx? - if i 



X' X 



{4 - l) £ , (46) 



<i - (^f) 2 fl4 W " O a 5 • (47) 

- 2 - ^ Qi ^" 4 (^ + ^-2) 2 ^ , (48) 



a xy = I — I a 



27 (^) « 4 (*? + *? - 2 



^ = 2 7 ^ « 4 K + , 2 - 2 , (49) 



2 



< = 2 7 (^) 2 « 4 (*? + 4 - 2) a , (50) 



where 
and 



Xi = x/a , ?/! = y/a , Z\ = z/a , (51) 

According to (|39|), what we shall need is the sum + 2 W^ 2w ) . To get it, we shall 

integrate over the volume: 

+ 2iy( 2 ^) = dx J a dy J C dz [dW^/dV + 2 tW^/^} = 

= a 2 c ^ darj ^ rf^i / x ^ ^ { ~ ^ ( Tr a ) 2 + a ^ + a L + 2cT i + + -( 53 ) 
Plugging of ( fi"5D - (|50[) in the above expression entails: 

W M + 2W^ = 2" 5 ^ (pfii) 2 ^ in ^ d Vl J\ d Zl { - i (x 2 - y 2 ) 2 + 



+ ~ (*? - i) 2 + ~ (</ 2 - i) 2 + (4 + v\ - 2) 2 + 

+ 2 7 (A + z 2 - 2) 2 + 2 7 (y 2 + z 2 - 2) 2 } « 

{pfllY 2~ 5 (63 7 + 20)(54) 



a 6 c 



fl 

where the numerics was performed by means of the Maple software. The product p 2 Q\ 
emerging in the above expression may be cast in a form wherein its temperature-dependence 
becomes manifest. According to (E7]), (EH), and (p9|), 



P 2 «i = M ^ 7 09 ^s) 2 sin 4 # (55) 

[1 + (c/a) 2 J 4 
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where k signifies the Boltzmann constant and (5 stands for the parameter of suprathermality: 

^ = ^ = PkT gas (56) 

The cosmic dust is thermal when (3=1. This way of writing the energies enables one to 
express these in terms of the temperature of the surrounding gas, T gas : 

W (w) + 2 = ^ ^ (63 7 + 20) - °" ((3 kT ga8 f sin* 8 (57) 

fi [1 + {c/a) 2 \ 



de/dt = - ° - ™> r _ „ (/3 fcT gas ) 2 sin 3 ^ (58) 



Together, (|38|), @ and (0) give: 

3 2 2~ 7 (63 7 + 20) a~ 2 c~ l 
liQJ [1 + (c/a) 2 ] 

According to (^) and (|27|) , 

/39 X 1 / 2 

J = {2 [3 k T 9as h) 1 ' 2 = (y /* * T gas p a 4 c) = 2 5 / 2 S" 1 / 2 a 2 c 1 / 2 p 1 / 2 {(3 fc T, as ) 1/2 

(59) 

Substitution of the latter in the former will give us the final expression for the alignment 
rate: 

5 - 5 ^\ 3/2 ,„ , m x3/2 sin 3 # 



a a 



M/dt = " U {P JQW> 3 2 ^ 63 (e/a) cot ° + 20 

(60) 

Simply from looking at this formula one can conclude that the major- inertia axis slows 
down its alignment for 9 approaching zero. This feature looks physically reasonable. 

On physical grounds, one may also expect that the alignment rate vanishes for h 
approaching unity (i.e., when the body lacks an axis of maximal inertia). Besides, one may 
expect the major axis of the grain to linger in the position 9 « 71 /2 , as if "hesitating" 
to whether to start aligning along or opposite its angular momentum. However, none of 
the latter two features seems to be instilled into (|60"D : it may seem from this formula that 
d9/dt is h— independent, and that the major axis leaves the initial position 9 = tt/2 with 
a finite angular velocity. To dispel these illusions, simply recall that our treatment is valid 
only for as long as the rotation and precession are fast motions compared to the alignment: 
9 < 00 . 

All these subtleties are anyway irrelevant when one merely wants to estimate the 
relaxation time, i.e., the time required for the maximal- inertia axis of an oblate cosmic-dust 
grain to be considerably shifted toward alignment with the angular momentum: 

'd9\ V 1 r s d9 



tl ~ l\dt) J ~ L/2 de/dt ' (61) 
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where 5 is introduced to avoid the divergence associated with the "slow finish". (One can 
take, for example, 5 = tt/8 .) A particular choice of S will bring into the expression for 
t\ some numerical factor of order unity. Since we want nothing more but a rough estimate 
for t\, we shall approximate t\ simply by the inverse d6/dt evaluated in the middle of the 
interval, at 9 = tt/4. Then 

U « <f [1 + (c/aff (£f U, ^ , Q P ^ -g^ (62) 

The quantity (c/a) 3 / 2 [1 + (c/a) 2 ] 4 [63 (c/a) 4 + 20] _1 is a steep function of parameter 
(c/a) . For c/a = 1/2 it equals to 4.3 x 1CT 2 , while for c/a = 1/10 it will be 
1.6 x lO- 3 . 

The above formula for the typical time of alignment is the main result of our article. 
Now we must think of the possible values for the material parameters involved. The 
values may depend both on the temperature T grain of the cosmic-dust grains, and on the 
frequency of the precession. To start with, temperature- and frequency-caused variations of 
the density p may be neglected, as they are to be small anyway. This way, we can use the 
(static) densities appropriate to the room temperature and pressure. 

Now, consider the isothermal shear modulus fi . The tables of physical quantities 
would provide its values for the room temperature and atmospheric pressure, and for 
quasistatic regimes solely. As for the possible frequency-related effects in materials (the 
so-called ultrasonic attenuation), these become noticeable only at frequencies higher than 
10 8 Hz (see section 17.7 in Nowick and Berry 1972). Another fortunate circumstance is 
that the pressure- dependence of the elastic moduli is known to be weak (Ahrens 1995). 
Besides, the elastic moduli of solids are known to be insensitive to temperature variations, 
as long as these variations are far enough from the melting point. The value of fi may 
increase by several percent when the temperature is drops from the room temperature to 10 
K. Dislocations don't affect the elastic moduli either. Solute elements have very little effect 
on moduli in quantities up to a few percent^. As for the role of the possible porosity, the 
elastic moduli scale as the square of the relative density. For porosities up to about 20 % , 
this is not of much relevance for our estimates^. 

According to (Ryan and Blevins 1987), at T ~ 20 K , the share modulus value 



Beyond that, one might assume that the moduli vary linearly with substitutional impurities (in which 
the atoms of the impurity replace those of the hosts). However hydrogen is not like that: it enters the 
interstices between the atoms of the host, and has marginal effect on modulus. 

5 We are deeply thankful to Michael Ashby and Michael Aziz, who consulted us on all these subtle topics. 
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for silicate \l ~ 10 10 Pa and density is p slllcate « 2500 kg m~ 3 . The temperature of 
subthermal rotating that according to Lazarian and Draine (1997) is important during 
crossovers is T rot ~ 20 K . 

Now, several words on the choice of values of the Q— factor will be in order. Purcell 
refers to (Krause 1973) wherein a review of acoustic dissipation in silicates is presented. 
According to (Krause 1973), for frequencies varying from 50 kHz to 27 GHz , and 
temperatures varying from 10 K to 50 K , the measured values of Q range between 
400 and 2000 . We would tend to believe that interstellar grains are almost certainly have 
plenty of cracks and defects. As known from seismology, for real silicate rocks the Q— factor 
is typically between 150 and 300 . Therefore we believe that Q( sl1 ) < 400. We are not 
familiar with the measurements of Q for carbonatious materials, but it is likely that Q for 
them will be lower than for silicates. 

Therefore, for the silicate grains with axis ratio 1:2 we get: 

( < 53 8 x 108 w (mT ' (63) 

while for graphite grain we assume that the axis ratio is 1:10 Q — 100 and obtain 

<■ * 7 x 106 w fe)" • (64) 

where the values of a are supposed to be in metres while t is in seconds (so that for 
example for a = 10~ 7 m , c/a = 1/2 and (3 = 1 the typical time will be 8 x 10 9 s ) . 

We have not attempted to calculate the dissipation rates for grains formed by loose 
aggregates of smaller particles. The inelastic relaxation within such grains may be orders of 
magnitude more efficient due to friction between parts of it ("effective viscosity"). It can 
well dominate the Barnett relaxation even for much larger grains. 



7. Dynamics of Prolate Cosmic-Dust Grains: Libration. 

At firts glance, the dynamics of a freely-spinning prolate body obeys the same 
principles as the dynamics of an oblate one: the axis of maximal inertia will tend to align 
itself parallel to the angular momentum. If we model a prolate body with a symmetric 
top, it will be once again convenient to choose it be a prism of dimensions 2a x 2a x 2c, 
though this time half-size c are larger than a , and therefore I3 = I2 > I\ . Then 
all our calculations formally remain in force, up to formula (|3§): since now the factor 



h — 1 = [1 — (c/a) ]/[l + (c/a) ] becomes negative, the right-hand side in (|36|) will change 
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its sign: 



dE rot J 2 

ro = — (h — 1) sin 9 cos 9 = — uj J sin 9 . (65) 
d9 J3 



Thereby formula will also acquire a "minus" sign in its right-hand side: 



-1 



d9 I dE rot \ dE rot , , . /n -i 

~dl ~ [ d9 I dt 



lo J sin 9y l W (66) 



Formulae ([J?p and (|39|) will remain unaltered. Eventually, by using ( pB| ) and (|3"9"D, we shall 
arrive to a formula that differs from ([51]) only by a sign, provided we keep the notation 9 
for the angle between J and the body-frame axis 3 (parallel to dimension 2c): 



/ a\ 3/2 ln lrri ,3/2 sin 3 # 



= " jrfmr (/3 ^" )3/2 325 2 "'' 5 (63 (c/a)4 cot2fl + 20 

(67) 

This looks like axis 3 tends to stand orthogonal to J, which seems to be so natural since 
axis 3 is now not the maximal-inertia but the minimal-inertia axis. 

Alas, all this nice extrapolation of the oblate-body-applicable approach to a prolate- 
body case is of no practical interest, bacause in reality an innnitesimally small deviation 
between the values of I2 and I3 leads to a considerably different type of wobble: the 
so-called libration (Synge & Griffith 1959). This phenomenon will be comprehensively 
discussed in our next article. 



8. Discussion 

8.1. Comparison with the Barnett relaxation 

Another important mechanism of internal relaxation, i.e. the Barnett relaxation, 
dissipates the energy via oscillating magnetization that arising from angular velocity 
precession in grain body coordinates. Lazarian & Draine (1997) provide an estimate 
t B ~ A/uj 2 , where A m 7.1 x 10 17 (a/10" 5 cm) s for an oblate grain with 2:1 axis ratio. 
For such grains the Barnett dissipation dominates for grains with a < 2 x 10 -6 m. However, 
this is not true for a larger axis ratio. Indeed, according to Lazarian & Draine (1997) for 
a/c 1 t B scales as (a/c) 6 , while we found above that ti scales as (c/a) 3//2 . Therefore fore 
grains with 4:1 axis ratio the inelastic relaxation dominates if a > 10~ 7 m. For grains less 
than 3 x 10~ 8 m the inelastic relaxation becomes more important. It also dominates for 
grains constituted by loose aggregates of smaller particles. The exact value of the Q-factor 
and therefore of relaxation rates would depend on the structure and properties of the 
aggregate. 
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For grains with axis ratio 2:1 and radii 10~ 7 m < a < 2 x 10~ 6 m the Barnett relaxation 
is dominant when grains rotate thermally. Suprathermally rotating grains stochastically 
undergo spin-ups and spin-downs and during crossovers rotate with subthermal velocities. 
Although short, in terms of grain alignment, crossovers, are the most important moments 
of grain dynamics. Grains are marginally susceptible to the randomization via gaseous 
bombardment when they rotate suprathermally (Purcell 1979). Our calculations of the 
inelastic relaxation efficiency show that the Barnett relaxation is the dominant process that 
determines internal dissipation during crossovers within grains with a > 1CT 7 m. Therefore 
results on paramagnetic and mechanical alignment obtained for such grains in Lazarian & 
Draine (1997) stay unaltered provided that Q factors of interstellar grains are as high as 
they were chosen in this paper. 

In this paper we assumed that grains are axially symmetric (and oblate). For grains of 
arbitrary shape accelerations will be higher as it tumbles. Therefore the inelastic relaxation 
is bound to increase and a relaxation time to decrease by a numerical factor that will 
depend on the grain shape. This subject is beyond the scope of the present paper. 

Our treatment of the inelastic relaxation ignored thermal fluctuations in grain 
material. In reality, for finite grain temperatures, 9 fluctuates and the thermal distribution 
proportional to the Boltzmann factor exp(—Eki n (6)/kT) is established as t — > oo. To 
describe the transient processes of alignment one can solve Fokker-Planck equation as it 
is done in Lazarian Sz Roberge (1997) in the case of the Barnett relaxation, but to use 
coefficients derived in Appendix E. 

8.2. Alignment and Internal dissipation 

Several mechanisms of cosmic-dust alignment are known. They constitute three major 
types: mechanical, paramagnetic and via radiative torques. All of them appeal to internal 
relaxation that enables the alignment of the angular momentum with the axis of the 
maximal moment of inertia, henceforth the axis of major inertia. The degree of achievable 
alignment depends on whether grains rotate thermally f] (the average rotational energy 
of a grain is of the order of the kinetic energy of the surrounding gas) or suprathermally 



6 In Landau and Lifshitz (1969, section 26) it was shown that the equilibrium position of an inelastic body 
corresponds to its maximal-inertia axis being parallel to the angular momentum. However in reality this 
statement remains true only up to thermal fluctuations: the rms of the angular deviation is A9 — \J IkT/J 2 , 
J and I being the angular momentum and the moment of inertia (Lazarian 1994, Lazarian and Roberge 
1997). 
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( spinning with energies much exceeding kT). 

We note that all the types of mechanisms below provide the alignment of the grain 
angular momentum in respect to the magnetic field. This happens because grains swiftly 
precess about magnetic field lines. This precession is called into being by the interaction of 
the grain's magnetic moment with the field. The said magnetic moment is generated by the 
Barnett effect and is thereby parallel to the angular velocity. For a typical cosmic-dust grain 
(~ 10~ 5 cm) in a typical interstellar magnetic field (5 /iG) the period of Larmor precession 
is less than a week (Purcell 1979), which is much less than the typical time of alignment. 

To relate the alignment of angular momentum to the polarization produced by 
grains, one has to know the alignment of grain axis, which is determined by the internal 
dissipation. However, the alignment of angular momentum, in its turn also depends on the 
alignment of grain axis. In the cases of thermally rotating grains, e.g. Davis & Greenstein 
(1951) and Gold (1952) alignment mechanisms it is sufficient to know that the rate of 
internal dissipation exceeds the rate of alignment. In cases of suprathermally rotating 
grains more precise estimates of the internal relaxation time are necessary. The latter 
include paramagnetic alignment of suprathermally rotating grains (Purcell 1979, Spitzer 
& McGlynn 1979, Lazarian & Draine 1979), crossover and cross sectional mechanical 
alignment (Lazarian 1995, Lazarian & Efroimsky 1996, Lazarian, Efroimsky & Ozik 1996) 
and the alignment via radiative torques (Draine & Weingartner 1996, 1997). Our finding 
that for some grains inelastic relaxation dominates is important and will be accounted in 
the quantitative studies of alignment elsewhere. 



8.3. Comparison with the earlier work 

Our treatment differs from that in (Purcell 1979) in a number of points. 

I. We have obtained different expressions for stresses. In particular, our stresses do 
obey the boundary conditions: they vanish on the grain boundary (while the stresses in 
(Purcell 1979) fail to do so). 

II. We have found energy dissipation at the double frequency. It is possible to show 
that the double-frequency dissipation provides a leading contribution []. The dissipation at 



7 For example, in formula ( |67| ) the term 63(c/et) 4 cot 2 9 in the sum in brackets is due to damping at the 
precession frequency uj , while the term 20 is due to damping at the double frequency. Evidently, even 
for c/a s 1/2 the double frequency provides the leading input, while for c/a w 1/10 its input becomes 
absolutely overwhelming 
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the double frequency was not considered in (Purcell 1979). 

As a result, our final expression for the same Q-factor predicts higher inelastic 
relaxation efficiency than the corresponding one in (Purcell 1979). To make a comparison 
we rewrite our formulae in a form similar to that in Purcell (1979). Then our expression 
(|60|) would read: 

- - I - " gP f . , 2l4 (0-4 (c/a) 4 cosH sine + O.l) (68) 
5 uQ 1 + (c/a) 2 v ' 

_ ir result ' w L 

In the appropriate expression presented in Purcell's article, the second term in brackets 
was missing (because Purcell missed the contribution from the dissipation at the double- 
frequency mode), while the first term is about 2.5 times larger (because Purcell's calculation 
of the stress tensor ignored the boundary conditions upon stresses). One may see that 
the our calculations reveal that the inelastic relaxation being much more efficient than 
predicted in Purcell (1979). For example, if we define a typical alignment time t« as the 
value of (dO/dt)^ 1 at 6 = tt/4 then, for silicate grains with c/a ~ 1/2 , "our" typical 
time derived from fl6~8"|) will be 23 times less than the appropriate estimate^] derived from 
Purcell's formula. For graphite grains with c/a ~ 1/10 , "our" typical time will be almost 
1200 times shorter than the estimate in Purcell (1979). 

III. In reality, the difference between our calculation and the one by Purcell will be 
even bigger: due to imperfectness of the grain material, we would choose for the Q— factor 
lower values than the one suggested by Purcell. (See the end of Section VI for discussion.) 



9. Conclusions 



The principal results of this paper are as follows: 



I. The stresses arising from tumbling of a grain deform it and change its moment 
of inertia, but this change lags behind grain angular velocity due to non-elastic effects. 
This process causes alignment of grain angular momentum with the grain axis of maximal 
moment of inertia. 

II. Deformations in tumbling grain happen both at the frequency of precession and at 
a double frequency. 



8 Comparing our results with that of Purcell, mind that the value of U obtained in (Purcell 1979) was 
for the degree of suprathermality (3 — 100 . 
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III. Inelastic relaxation dominates the Barnett relaxation for (1) grains with large axis 
ratio, (2) grains with low Q factor, i.e. fractal grains produced by coagulation, (3) grains 
with a > 2 /im, (4) suprathermally rotating grains. 
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A. Intermediate system of reference 

As pointed in Section III, our expression (|T6| ) for the acceleration of a point inside the 
spinning body considerably differs from an appropriate expression presented in the article 
(Purcell 1979). To be confident in our result, we shall now reproduce it by means of a 
two-step calculation. 

To do so, we shall introduce an "intermediate" coordinate system (X', Y', Z'), with 
basis vectors denoted as e^', ey/, and e z >. This system is obtained from (X, Y, Z) by 
means of consequent rotations of axis X by the Euler angles ip and 9. 

In terms of the said angles, 

Qi = <p sin 9 sin ip + 9 cos ip , (Al) 

Q2 — <fi sin 9 cos ip — 9 sin ip , (A2) 

Q 3 = (j> cos 9 + ip , (A3) 

and 

n ± 2 = ttj +nj = 2 sin 2 # + 9 . (A4) 

The position, velocity and acceleration of a point inside the grain, relative to the 
"intermediate" frame, will be denoted by r', v' and a'. Evidently, 

r" = r' = r , e 3 = e z , . (A5) 

Now, one will be able to calculate the acceleration a as 

a = a' + (e* <j>) x r' + 2 (e, <p) x v' + 2 e z x (e z x r') (A6) 

and plug, instead of a' , the expression 

a' = a" + (e 3 x r" + 2 (e 3 ^) x v" + V> 2 e 3 x (e 3 x r") (A7) 

This will once again lead to (|16|) . 



B. Diffusion coefficients 

In the absence of gaseous bombardment grain angular momentum stays constant. 
The alignment of J in grain axes is being determined by internal relaxation and that tend 
to decrease 9 and thermal fluctuations that randomize 9. To quantify these processes 
Fokker-Planck equation may be used. 
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In the spherical coordinate system J, 6, , where the polar axis is parallel to the 
principle axis of maximal inertia, the Fokker-Planck equation will read: 

f = -V.S, (Bl) 

Here / = /(J) is the joint distribution for J, 9 and , while S is the probability 
current: 

S = Af - -V -(Df) . (B2) 

where 



is the mean torque, and 




A - 77 > (B3) 



/AJAJ\ 

s (— > • < B4 > 



where (..) denote ensemble averaging, is the diffusion tensor. Generally speaking, A and 
B include the cumulative effects of all the processes that change J in the body frame. 
Still, in what follows we shall consider the inelastic relaxation as if it were the sole factor 
contesting the gas damping. All the other mechanisms of orientation will be ignored. 

We have written down the transport equation, and introduced the entities A and 
B, in agreement with (Risken 1984). The difference between our definition of A and the 
definition presented in (Landau and Lifshitz 1981) stems from the fact that in (Risken 
1984) the change of momentum of the grain is denoted as p + q , whereas in (Landau and 
Lifshitz 1981, formula (21.1) and thereabout) it is denoted by p — q . For this reason A 
in formula (21.5) for S, in (Landau and Lifshitz 1981), appears with a negative sign. 

The inelastic relaxation leaves J unaltered, in the inertial frame. In the body frame, 
which we are using here, the direction of J will vary, but its absolute value will be conserved. 
For this reason the components of A and B, having J in their subscripts, will vanish. As 
for ^-dependence, it would be nonexistent should the grain be an oblate ellipsoid rather 
than a prism. Still, in our approximation we may neglect the ^-dependence: it is, in fact, 
not that relevant for our purposes whether the grain a circular or a square cross section. 
That is, for the purpose of our estimates (which are anyway not exhaustingly exact) we 
feel free to consider the grain either square, when we need to estimate t a and d9/dt , 
or circular, when we need to simplify the Fokker-Plank equation. In other words, let us 
assume that in our study of the relaxation over 6 the ^-dependence is to be averaged out. 
Hence all the components of A and B with <p in their subscripts will be thrown out either. 
Now the only remaining components will be Ag ^ and Bgp where the extra superscript 
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is introduced to emphasize that only the effect of inelastic relaxation (versus the gas 
damping) is taken account of. 

The expression for A g l ^ is straightforward from our expression for d6/dt : 

AP = JdO/dt (B5) 

while BgJ' is to be found from the principle of detailed balance: in case the distribution 
function /(J) is the thermodynamical-equilibrium one ( / = Jte ), the rate of every 
microscopic process equals that of its time- reversed counterpart, and thence the probability 
current must vanish at every point in the phase space. In particular, 

= S e = sm6 f TE 4 r) - JL A { sm9 fTEBfp} (B6) 

where the thermodynamical-equilibrium distribution function is a Boltzmann one[]: 

f TE (J) = C exp{-E rot (e,J)/kT g ^ in } , (B7) 

C being a normalization constant, T grain being the grain temperature, and E rot (9, J) being 
the grain rotational energy (expressed by equation fl37|)). The solution to this first-order 
differential equation is: 



ir 2J 2 exp(esin 2 g) f . 2 , ~ 2 - . 2 x 2 J 2 A h 

B ee = , ■ a / sin y (A cos y + Bsm y)exp{-£x )dx — exp(- -4 

U sine Jtt/2 '• ' U a — 1 



2 ^ 

B8) 



where 



and 



, 2^63(c/ay 

A ~ 6WT20 ' (B9) 

B = -^^P (BIO) 

63(c/a) 4 + 20 V ; 



9 In the case of a grain being in an equilibrium, we formally treat it as a small system placed in a bath 
of temperature T gra j n . This is certainly fair, since the number of rotational degrees of freedom of the grain 
(= 3) is much less than the amount of its vibrational degrees of freedom (which is of the same order as the 
number of atoms in the grain) . 
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The solution ( |B8|) is obtained assuming that -B^(vr /2,£) = — J 2 
is smooth at 7r/2 (compare with Lazarian & Roberge 1997): 



A/ti, to insure that B\ 



'00 



■ir 



lim 



rf 2 ^(7T/2,e) 



< oo 



(B12) 



d6 2 



The coefficients and B l e r e can be used to describe internal alignment in the presence 
of gaseous bombardment and inelastic relaxation as it has been done in Lazarian & Roberge 
(1995) for the case of Barnett relaxation. 
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Fig. 1. — The coordinate system (e x , e y , e z ) is associated with the inertial frame, so that the 
(conserved) angular momentum J is aimed along e z . The coordinate system (ei,e 2 ,e 3 ) is 
associated eith the three principal axes of inertia of the body, so that e 3 points along the 
axis of major inertia. 



